The LU-LC conjecture is false 



alent stabilizer states are also local Clifford equivalent. If this conjecture were true, the local 
equivalence of stabilizer states would be extremely easy to characterize. Unfortunately 
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The LU-LC conjecture is an important open problem concerning the structure of entan- 
glement of states described in the stabilizer formalism. It states that two local unitary equiv- 

O' 

<n : 

however, based on the recent progress made by Gross and Van den Nest, we find that the 
conjecture is false. 

o : 

m . 
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43 ■ 

The stabilizer formalism is a group-theoretic framework originally devised to systematically 
analyze various quantum error-correcting codes [6] and has found applications in other areas of 
quantum information processing, such as the one-way computation model A Ell] and the quan- 
i— i , turn sharing of classical secrets HQ. 

I One of the key features of stabilizer states defined by the stabilizer formalism is the presence of 

high-degree multipartite entanglement in them. A natural approach to studying the properties of 
\Q ', the entanglement in stabilizer states is to investigate their local equivalences. Mainly, three types 

of local equivalences have been studied with respect to stochastic local operations and classical 
, communication (SLOCC), local unitary (LU), and local Clifford operations (LC). In the following, 

we will call two states SLOCC (LU, LC) equivalent if they can be transformed to each other by 
O ■ means of SLOCC (LU, LC) operations. 

. Local Clifford operations are local unitary operations that map the Pauli group to itself under 

^ ■ conjugation. Therefore, the action of a local Clifford operation on stabilizer state can be repre- 

sented as a simple way of translating the stabilizers. Because of this close relationship between 
the stabilizer formalism and the Clifford group, the LC equivalence of stabilizer states has been 
studied thoroughly in the literature. For example, a polynomial time algorithm has been found to 



decide whether two stabilizer states are LC equivalent [14] and the action of local Clifford group 



on graph states (an important subset of stabilizer states) has been translated into elementary graph 



transformations characterized by a single rule [15]. Specifically, the classification of LC equivalent 



stabilizer states has been performed systematically up to 12 qubits ESQ. 

Stabilizer states have extremely symmetric structures which put strong restrictions on the local 
unitary that can map one stabilizer state to the other. In fact, it has been conjectured for several 
years that any two LU equivalent stabilizer states must also be LC equivalent (the LU-LC conjec- 
ture, listed as the 28th open problem in quantum information [12]). If this conjecture were true, all 



three local equivalence of stabilizer states would be the same, as it has already been proven that 



2 



two stabilizer states are SLOCC equivalent if and only if they are LU equivalent [13]. Moreover, 
deciding whether two states are locally equivalent would be efficient and the local equivalence 
of stabilizer states could be described as purely graph theoretic terms. The conjecture has been 



18T) and is further supported by the 



proved for large subclasses of stabilizer states in Refs. 11/ 
results obtained in Ref. 11161] - 

In this paper, however, we will present a counterexample to the LU-LC conjecture which dis- 
proves the conjecture in general. Our result is build upon the recent progress that transforms the 
conjecture to a simpler problem [8]. In addition to our heuristics for the construction of coun- 
terexamples, we also present some of the special cases in which LU and LC equivalence are the 
same. 

The organization of the this paper is as follows. In the next section, we introduce the basics 
of stabilizer formalism and some notations and results used in this paper. Sec. [Ill] develops new 
representation of the problem in the first part and then considers some special cases of the LU-LC 
conjecture. The procedure of generating counterexamples and one explicit example are given in 
Sec.|lVl We conclude in Sec. [V] and discuss some possible future work on this problem. 

II. NOTATIONS 



Stabilizer states are quantum states described by a set of commuting operators. This idea 
of representing states with operators has been proved to be extremely useful in the theory of 
quantum information. Let {/, X, Y, Z} be Pauli matrices, 
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Q be the Pauli group generated by them, and J n be I® n . Mathematically, a stabilizer of n qubits 
is an Abelian subgroup of Q® n that does not include —I n . When the subgroup has cardinality 
exactly 2 n , there will be a unique quantum state determined by it as the simultaneous fix point of 
all the operators in the subgroup. For example, the EPR pair (|00) + |ll))/\/2 is stabilized by the 
group generated by {X (g) X, Z (g) Z}. For any graph G, the corresponding graph state is a special 
stabilizer state with 



X v Z u , 

u£N(v) 



v £ V(G) 



as its stabilizer, where N(v) is the neighbor set of vertex v. It is known that any stabilizer state is 
LC equivalent to some graph state (la] . Therefore, the LU-LC conjecture for stabilizer states and 
graph states are the same problem. For more details on the power of stabilizer formalism and 
graph states, the readers are referred to Refs. 

There are special structures in the amplitudes of stabilizer states expanded in the computa- 
tional basis. It is proved that any stabilizer state can be written as 

1 



T\ ^ 

1 x&T 



(1) 
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where T is an affine space of F?>, l(x) is linear in Xj with addition modulo 2, and q(x) is a quadratic 
function of x/s. Conversely, any state having the above form is a stabilizer state [3]. 

The main topic in this paper is to study local equivalences of stabilizer states, especially the 
LU equivalence and LC equivalence. Clifford operators are 2 by 2 unitary operators U such that 
UQU^ = Q. Up to a global phase, the Clifford operators form a finite subgroup of 17(2). We say 
two stabilizers l^o) and \ip\) are LC equivalent if there exists Clifford operators Uj such that 



3=1 

In the following, we define several problems, each of which is given an abbreviated name for 
later references. The first one is the original LU-LC conjecture itself, then a restricted case of it, 
DLU-LC. The third is the DLU-DLC problem and the last is quadratic form phase problem (QFP). 
One can see the close relationship between QFP and the other three problems from Eq. (d). The 
previously known relation of these problem is as follows. DLU-LC is shown to be equivalent 
to the LU-LC conjecture uM- DLU-DLC and QFP are equivalent and they imply the LU-LC 
conjecture according to Ref. [8]. All of these four statements are false as indicated by our coun- 
terexample given in Sec.HVl 

Problem 1 (LU-LC). Every two LU equivalent stabilizer states are also LC equivalent. 

Problem 2 (DLU-LC). Every two stabilizer states that can be mapped onto each other by means of a 
diagonal local unitary, are LC equivalent. 

Problem 3 (DLU-DLC). If two stabilizer states can be mapped onto each other by means of a diagonal 
local unitary, then also by a diagonal local Clifford operation. 

Problem 4 (QFP 18]). Let Sbea linear subspace o/F??, and Q : F% — > F2 be a quadratic function. If there 
exists complex phases {cj} such that 

n 

(_!)<?(*) = JJ c *i ( f or every x e S , (2) 
3 

then the phases can be chosen from {±1, ±i}. 



III. SIMPLIFICATION OF THE QFP PROBLEM 



Let us start by simplifying the QFP problem. We will give a linearized version of QFP in 
the following. It's easy to obtain linear representations of QFP by, for example, enumerating all 
possible x G S and resulting in a system of linear congruence equations. But this method does 
not give us much information on how to tackle the problem. What we present in the following is 
a more symmetric linear representation that leads to the construction of counterexamples of QFP. 
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A. Symmetric Linear Representation 

There are two main mathematical objects in the problem, the quadratic function Q(x) and the 
subspace S. We will find a representation of Q(x) first and then analyze the structure of S. Based 
on this analysis, some special cases of QFP will be proved. We will also obtain an approximate 
restatement of the general QFP problem that leads to the construction of counterexamples. 

First, noticing that the linear terms in Q(x) can be moved to the right hand side of Eq. (O 
without changing the problem, we can consider quadratic forms only. There is a trivial one-to- 
one correspondence between quadratic forms over F2 and simple graphs. For each Q(x), one can 
identify it with a graph Q with vertex set V = {1, 2, . . . , n} and edge set 

E(Q) ={(i,j) I XiXj is a term in Q(x)}. 

Each subset A of V defines a subgraph Q\j± whose vertex set and edge set are A and E(Q) n A 2 
respectively. With the above definitions, one can verify that 



Qip)= E{Q\j) 



(3) 



where I x is the indicator set {j \ Xj = 1}. 

Next, to analyze the structure of S, we choose a basis j^ 1 , £ 2 , . . . , £ d } of S. For each x <E S, we 
can represent it as the linear combination 

d 
k=l 

Let h £ be the vector whose k-th coordinate is and we label the above x with it as x h . The 
size of S is therefore D = 2 d . 

We introduce a new concept called pattern of a position in S. The pattern for position j is de- 
fined to be the vector meF^ whose k-th coordinate equals to See Fig.[T]for a illustration of 
this definition. Patterns are basis dependant but the choice of basis does not change the following 
analysis. 

Let A m be the set of positions having pattern m, that is, 

A m = {j I m k = (* for all k}. 

The collection {A m } is obviously a partition of the vertex set V of graph Q. For any j G A m , one 
can calculate the j-th coordinate of x h as 



x j = Yl hk $ = Yl hkmk d = m ^ ' ( 4 ) 



k=l k=l 



where the summation is taken over F2, or modulo 2. That is, the value of x h at position j is 
determined solely by the pattern of the position. We can define variable xr m i to be a representative 
variable of those having pattern m. 
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FIG. 1: Illustration of patterns 

Consider a special case of QFP where A m is nonempty for all patterns mGFj- {0}. We will 
prove that in this case QFP is true and the details of the calculation will be useful even when we 
deal with the general QFP problem. 

As all A m 's are nonempty, we can rewrite Eq. (O in the QFP problem by replacing the variables 
with their representatives in the right hand side, 



= n c; 



for every x G S, 



where Cr TO i is the product of all Cj for j £ A m . For a specific x h G S, this means that 



From Eq. (O and the fact that 



we have 



m:(m,/i)=l 



(5) 



(6) 



E 

m<m':(m,/i)=l,(m',/i}=l 



-En 



q(**) = E l^ccuj 

m:(m,h)=l 

= Y^{m,h) E{Q\ Am ) + ^ {m,h)(m,h)E mm <, 

tti m<m' 

where -E mm ' counts the number of edges between vertex sets A m and A m > in graph Q. The first 
part of the above summation can be omitted as it can be absorbed into the right hand side of 
Eq. (O by changing the values of Cr m i appropriately. Equation ^ is therefore reduced into the 
following form, 



^2 (rn,h)(m',h)E mm/ 

(-1) m < m> = n c M h) - 



(7) 
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Let rr m i be the real number satisfying 

C M =?M. (8) 
By taking logarithm on both side of Eq. 10, one gets 

^2 h ) r [m] = 2 ^2 ( m i h)(m', h)E mm , (mod 4). 

mj^O m<m' 

This equation holds for all h G F^, ant it actually specifies a system of D — 1 equations where 
L> = 2 d . By defining D — 1 by D — 1 matrix G = and matrix T with element 

T itm = (i,j)(i,k), i,j,ke¥i-{0},j < k, (9) 

we can write the system of equations succinctly as 

Gr = 2Te (mod 4) (10) 

where r and e are vectors consisting of rr m i and E mm > respectively. 

Next, we calculate G^ 1 and the product 2G _1 T based on the following lemma. 

Lemma 1. The following identities hold with the indexes traversing — {0}. 

Y,{i,j) = D/2, (11) 

3 



, D/2 i = k 

52{i,j)m = { (12) 

DJA otherwise 



For k / I, we have 



D/4 i = kori = l 

i = k®l (13) 

£)/8 otherwise 



We prove the last identity only others are simpler and can be dealt with similarly. When i = fc 
or i = I the result follows from the second identity. In the case of i = k © I, there does not exist j 
so that (j, i), (j, k), (j, I) are all 1. What remains to show is when i ^ k, i ^ I, k / I and i ^ k®l. 
Consider a 3 by d matrix F with i,k,l as its rows, then it is a matrix of rank 3 under the above 
conditions, j contributes 1 to the sum when it is a solution of Fj = (1, 1, 1) T . The number of 
solutions is obviously D/8. 

With the above identities, it's easy to check that G^ 1 is given by 

G- 1 = -|(2(z,i)-l), (14) 
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which is obtained by replacing the zero elements in G with —1 and multiplying a normalization 
constant. More importantly the product 2G _1 T is still an integral matrix with 1 at {i, (J, k)) where 
i = j or i = k, — 1 where i = j © k and O's elsewhere. This immediately proves the QFP problem 
in the special case where A m 's are all nonempty as we can always choose r to be 2G~ 1 Te, an 
integral vector. The proof is completed without even referring to the condition that there exists a 
real solution of f for Eq. <flQ]) as promised by the QFP problem. 

However, the general case QFP problem is much more complicated and, in fact, false some- 
times. When A m is empty for some pattern m, we cannot introduce representative variable acr m i 
as there is no corresponding variables. But in order to make use of the above formalism, we keep 
the variables xr TO i for all those missing patterns m and, at the same time, fix the corresponding 
C[ m ]'s to be 1, or fix r[ m ]'s to be equivalently. Hence, the QFP problem is approximately reduced 
to the following statement in terms of linear congruence equations (LCE). We have intentionally 
simplified the problem by omitting the indication of places forced to be in f. This gives us a 
stronger statement than QFP, that is, the LCE problem implies QFP and any counterexample of 
QFP also invalidates LCE. Compared to the trivial linearization of QFP, we have actually used 
more variables to make the representation more symmetric. 

Problem 5 (LCE). If the following equation (Eq. (fT0]> j of f has a real solution, 

Gr = 2Te (mod 4) 

then it also has an integral solution that preserves all zero entries in the real solution. 



B. The Low-Rank Case 



In this subsection, we discuss the QFP problem with low-rank subspace S. Namely, we want 
to show that for d < 5, there is no counterexamples of QFP. As QFP is a sufficient condition for 
LU-LC, it follows that LU-LC conjecture holds for d < 5. 

We need only to prove that there is no counterexamples for LCE in this case. First, any solution 
f of Eq. (flOl) can be written as 

r = 2G~ 1 Te + 4G- 1 s, 
for some integral vector s, and vice versa. Substitute Eq. ((141) into the second part, 

f=2G^Te + ^Gs-^l, 

where 1 is the vector whose entries are all 1 and a is the summation of all entries in vector s. 
When d < 3 all parts in the summation are integral. If d = 4, only the last term can be half 
integral. However, when this happens, no entry in r is integral, and the LCE problem holds by 
choosing s = 0. 

Consider now the case of d = 5. Let [x] be the largest integer not exceeding x, and {x} be 
x — [x]. We claim that if f is a solution of Eq. ((TO) , the truncation [r\ is also a valid solution. It 
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suffices to prove that, for any integral s, 

G [4G _1 s] = (mod 4). (15) 

Employing Eq. ((14)) and a detailed discussion on the four possible cases of a (mod 4), we can 
reduce the above equation to 




= (mod 4), 



which in turn follows easily from the second identity of Lemma [T]by noticing that 




is a column of G. 

We have shown that it is impossible to construct counterexamples of the QFP problem when 
the rank of subspace S is less than or equal to 5. Naturally, the next step is to investigate the case 
of d = 6, which will lead to counterexamples. 



IV. RANDOM GENERATION OF COUNTEREXAMPLES 



Notice that the "zero-entry preserving" requirement is essential in LCE. Without this addi- 
tional requirement, LCE will be true. Therefore, we try greedily to find a real solution of r with as 
many zero entries as possible. As the number of nonzero entries will be the size of our counterex- 
ample, this approach also tries to minimize the size of counterexamples. 

Recall the key equation of the problem 

r = 2G- 1 Te + AG~ 1 s. 

As the first part of the right hand side is integral, we will at least need to find an integral vector 
s so that 4G _1 s has many integral entries in order to have a r with many O's. Conversely, once 
4G _1 s is chosen, we can try to cancel most of its integral entries by choosing e. And after s and 
e are fixed, we can calculate r and verify whether there is no integral solution that preserves zero 
entries in it. It is a backward approach that chooses in order, s, e, and r. More details on this 
approach are given in the following. 

We first choose s"o at random. But it generally does not give us many integral entries in 4G~ 1 s~o. 
We can employ the rounding technique once again previously used in the case of d = 5. For d = 6, 
the following equation similar to Eq. (fl5]> can be proved 

G[8G _1 sS]=0 (mod 8). 

Therefore, there exists an integral s% such that 

^ [8G _1 a6] = 4G-Vi (mod 4). 
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There will be generally more integers in AG s\ than in 4G , ~ 1 s~o and we can use si instead of sq. 

In the next step, vector e is chosen to cancel the integral entries in AG~ l s{. Thanks to the simple 
form of 2G~ l T, it is flexible to find e that fulfills the requirement and, in addition, satisfies that 
E mm i = when pattern m corresponds to an integral entry in 4G , ~ 1 «s~!. 

Vector e gives most of the useful information on the quadratic form Q(x), and we construct 
subspace S by deleting columns of G corresponding to zero entry positions in r. This is how a 
candidate counterexample is generated. 

Finally, we need to verify whether the Q(x) and S obtained here constitute a valid counterex- 
ample of QFP, that is, to verify that there do not exist Cj G {±1, ±i} satisfying Eq. (|2]). This can be 
done by a modified Gaussian elimination method described below. 

Let A be the matrices contains all vectors in S as its rows. The problem of whether a can all be 
chosen from {±1, ±i} is equivalent to whether an equation of the following form has an integral 
solution 



The main difficulty in performing Gaussian elimination is that in congruence equations we 
cannot perform division to normalize the pivot. Fortunately, for equations of modular 4, division 
operations can be avoided. When performing the elimination procedure, try first to find a row 
with odd pivot and if the pivot is congruent to 3, multiply 3 on it. Otherwise, if we cannot find 
a row having odd pivot, the division is also unnecessary as one can do elimination simply by 
subtraction. Using this modified elimination method, we can derive from Eq. ((16)) a contradicting 
equation as 



under the assumption of the existence of an integral solution. 

The above procedure to randomly generate different counterexamples is implemented and is 
available at http://arxiv.org/e-print/0709.1266 as a gzipped tar (.tar.gz) file. Here is one of the 
counterexamples it generates. 

In the example, n = 27 and the subspace S of F 27 is given by a set of basis as 

= 100010001010101010100011110 
£ 2 = 101010111001100000001010101 
f = 011001100111100111100110011 

( 18 ) 

£ 4 = 000111100000011001100001111 
f = 000000011111111000011111111 
£ 6 = 000000000000000111111111111. 

One can see that every two columns (patterns) are different and the columns are organized in an 
increasing order. S has a rank of 6 and consists 64 elements. 
The quadratic form Q(x) is the summation of 11 terms: 

X\X2 + X1X3 + XlX 8 + X2X4 + X 2 X 8 + X 2 Xiq + X3X4 + X 3 Xs + X3X16 + X4X8 + XgXlQ. (19) 



Ax = b (mod 4). 



(16) 



= 2 (mod 4), 



(17) 
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We note that this is not directly generated by the random procedure. In fact, we have find a Q 
with smaller number of terms that has the same value Q(x) as the generated quadratic form for 
all x G S. Originally the quadratic form contains more than 60 terms. This simplification is 
suggested by one of the referees of QIP 2008. 

The QFP problem defined by the above Q and S has a solution where the phases are powers 
of e™/ 4 , but not in {±1, ±i}. We give the exponents sequentially: 

e = [3, 5, 7, 5, 1, 3, 5, 7, 1, 3, 5, 5, 3, 7, 3, 3, 7, 1, 7, 3, 1, 5, 5, 5, 3, 5, 3]. (20) 

Our final aim is to disprove the LU-LC conjecture and what we have already shown is that 
QFP is false generally. As proved in Ref. [8], QFP implies LU-LC. Therefore, QFP is false does 
not logically guarantee that LU-LC is false. Fortunately, however, the above counterexample of 
QFP can be transformed to a counterexample of LU-LC. Although not proved, this transformation 
seems to work all the time. 

Here is the way we perform the transformation. Define two states 



|*> = E 



iq,5) = E(-i) q(x)i ' 



x&S 



From Eq. (D, we know that they are stabilizer states. The corresponding QFP problem of the same 
Q and S indicates that they are LU equivalent, in fact, even DLU equivalent. We need to show 
that they are not LC equivalent. This can be done efficiently by using the LC equivalence decision 



algorithm |14I1 for graph states after efficiently finding graph states \Gs) and |Gq,s) which are LC 
equivalent to \S) and \Q,S) respectively. If the algorithm tells that \Gs) and \Gq : s) are not LC 
equivalent, nor are the states \S) and \Q,S) (see Fig. [2]). The four states in the diagram are in the 
same equivalence class under LU criterion. But when LC equivalence are considered, the upper 
two and the lower two belong to different classes. 



\Q,S) = ^= y £(-l) < *M\x) ^ \G. 
V l 5 l xes 



,si 

LU I NOT LC 



I 

\S) = -^=Y\x) ^ \G S ) 



FIG. 2: State Transformation 



A better solution to show the non-LC equivalence of \Q,S) and \S) would be a proof that 
DLU and LC imply DLC. This is open generally, but it is pointed out by B. Zeng that when local 
unitaries are all non-Clifford, just as in the above case, a proof can be found using results on 
minimal support of stabilizers developed in Ref. [17]. 

The two graph states \ Gs) and \Gq,s) have corresponding graphs shown in Fig. 13 The graphs 
with and without the dotted edge are graph Gq : s and Gs respectively. Interestingly two LU 



11 



equivalent but not LC equivalent graph states can differ only in on edge. Note that the simplifi- 
cation of Q(x) discussed above does not simplify or change the form of the four stabilizer states 
in Fig.lU 



In summary, we have shown that the LU-LC conjecture is false by giving explicit counterex- 
amples of it. It is also clear that in order to disprove the conjecture, we have to consider stabilizer 
states with the rank of the support no less than 6. This result leads us to rethink about the local 
equivalence problem of stabilizers as most of the previous work focus on proving the conjecture 
to be true. 

The random procedure generates counterexamples of size 27 and 35. Though not proved, we 
believe that 27 is the smallest possible size of counterexamples of LU-LC. Although the LU-LC 
conjecture is now disproved, we still know little about the relation of LU and LC equivalence and 
do not have an explicit understanding of why or when LU equivalence differs from LC equiva- 
lence. When d = 6, the random generation procedure can find a counterexample in seconds, but 
it is in fact not an efficient algorithm when d is large or even when d = 7. That fact is we have 
never obtained a valid counterexample of d = 7 using the random procedure. Fortunately, larger 
scale counterexamples, including those of d = 7, have been found [7] motivated by the randomly 
generated counterexamples. 

As LU and LC equivalences are now known to be different. It is also challenging to ask whether 
there is an efficient algorithm deciding LU equivalence for stabilizers or whether there is a graph 
theoretical interpretation of LU equivalent graph states. This seems to be difficult as local unitary 
operations are much less linked up with the stabilizer formalism. 




FIG. 3: Corresponding Graph States 



V. CONCLUSION 
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